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1. Introduction 


Non-leptonic weak decays of hyperons have been of interest for many years. Especially 
mesonic decays, such as A —>■ Ntt, are studied in order to reveal the properties of the low 
energy weak interactions among quarks. Experimental data for such decays indicate a strong 
AJ = 1/2 enhancement compared to AJ = 3/2 component. This AJ = 1/2 enhancement, 
known as the AJ = 1/2 rule, is not expected naively in the standard model of the weak 
interaction, and therefore its origin should be attributed for corrections to the weak vertex 
due to the strong interaction. Part of the strong corrections can be estimated by using the 
renormalization group improved perturbation theory of QCD[^^^,0, while contributions of 
the low-energy hadronic interactions are not quantitatively understood. 

In order to study the corrections to the weak interaction due to the low-energy hadronic 
interaction, it may be useful to look into a new type of weak processes, such as NN ^ NN, 
and AN —>■ NN, i.e., the two-body weak scattering processes. It is known that the AN —>■ NN 
transition plays a dominant role in the nonmesonic decays of hypernuclei, whose data have 
been accumulated in recent hypernuclear experiments. Therefore it seems timely to study the 
AN —>■ NN weak transition from the standard theory point of view. 

The purpose of this paper is to study the roles of quark structure in the two-body AN —>• 
NN weak transition, and to construct the induced transition potential Because the 

AN NN decay has a large momentum transfer of approximately 420 MeV/c (assuming 
the relative momentum of the initial A and N is zero), the short distance dynamics of two 
baryons must be signihcant. We here propose that the AN NN transition at short distance 
is described by a direct quark mechanism, where a contact four-quark interaction between the 
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Figure 1: The diagram for the one pion exchange mechanism 

constituent quarks of baryons causes the transition without exchanging mesons. The four- 
quark vertex is taken from the low energy effective weak Hamiltonian in the standard theory. 
It contains both the AJ = 1/2 and AJ = 3/2 components. We will see that no AJ = 1/2 
enhancement is seen in the two-baryon process, AN —>■ NN, and thus the AJ = 3/2 component 
of the weak interaction may be observed. 

We will compare the direct quark transition potential with the transition potentials based 
on the meson exchange mechanism (Fig. p, such as n, K, p exchanges 0000. There the 
meson-baryon-baryon (such as tt —A —iV) weak vertex is determined phenomenologically so as 
to describe the free hyperon decays, and therefore satisfies the AJ = 1/2 rule automatically. 
Recent experimental data, however, have revealed some difficulties in the meson exchange 
picture. For instance, the predicted n-p ratio is much smaller than the experimental data for 
light hypernuclei. 

This paper is organized as follows. In section 2, we present the effective weak Hamiltonian 
derived in the standard theory with perturbation in QCD. In section 3, the weak Hamiltonian 
is applied to the ANtoNN transition and the direct quark induced transition potential V{k, k') 
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is calculated. In section 4, we present the explicit form of the transition potential given in the 
momentum space. In sections 5 and 6, the transition potential is applied to the decay of light 
hypernuclei and the results are compared with available experimental data. Discussions and 
conclusions are given in section 7. 


2. Effective 4-quark weak interaction 


The standard model describes hadronic weak interactions by exchanges of the weak gauge 
bosons between quarks. Because the gauge bosons are very heavy, low energy phenomena can 
effectively be represented by a Hamiltonian composed of four-quark vertices, which contain the 
QCD corrections on the pure weak vertex. Such an effective Hamiltonian has been studied by 
several authors[0,|],^||. It can be computed by evaluating the perturbative QCD corrections, 
using the operator product expansion and the renormalization group equation for the Wilson 
coefficients. 

In the standard model, strangeness changing weak decay is described by the vertex 


Hw{x) 


9 

2v/2 


j;{x)w;{x) + H.c ., 


( 1 ) 


where is the charged-W-boson held and is the hadronic charged weak current. Note 
that the standard theory contains no havor changing neutral current. The effective Hamilto¬ 
nian for AS = ±1 nonleptonic processes is dehned by 

(I //fjf ‘ \) = -\jd‘x(\THw(x)Hw{0) |> (2) 


where RHS is the weak transition matrix element between low-momentum hadron states 
composed of light quarks and differing in strangeness by one. We separate the mass-scale 
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dependent coefficients and the four-quark operators using the operator-product-expansion. At 
the mass scale /i = Mw the QCD running coupling constant as is so small that the coefficients 
can be expanded perturbatively in ag- Paschos et al. 0 takes the following Hamiltonian at 


/i — Mw, 


H. 


AS=1 

e// 


fj,=Mw -^2 


[^u{SaUa)v-A{updf:i)v-A + ^c{SaCa)v-Aicpd0)v-A] + iPpen3(t-quark) (3) 


where a and (3 stand for color indices of the quark held. The hrst term is a pure weak 
interaction at low-momentum transfer with = VqdVq^ where matrix V is the Cabibbo- 
Kobayashi-Maskawa matrix. This term contains AI = 3/2 component as well as AI = 1/2 
component. The second term, TTpeng(t-Quark), represents the hrst order QCD correction that 
is produced by a so-called penguin diagram with the top quark in the intermediate line. This 
term is needed because the top quark is heavier than the W-boson. Starting from eq.(|^) 
at p = Mw the ehective Hamiltonian at low mass scale is computed with the help of the 
renormalization group technique. The one loop QCD corrections are taken into account. 
Operator mixing takes place and enhances the AI = 1/2 component, while the AI = 3/2 part 
is suppressed. 

The perturbation theory, of course, cannot be extended down to the low energy region, 
where non-perturbative ehects of QCD may modify the weak vertex as well. Here we employ 
the picture proposed by Bardeen et a/.[]^, i.e., we assume that the perturbative correction is 
applied down to /i^ ~ /Iq, where a{fil) = 1 The effective four-quark Hamiltonians is calculated 
at Pq, and is applied to the quark model calculation (or any other low energy theory, such as 
the chiral effective theory). 


5 






We use the following Hamiltonian in the present calculation : 


where 


and 




^ 6 

i: 

V ^ r=l,r/4 




Ki 

K 2 

7^3 

K, 

7^6 

-0.284 

0.009 

0.026 

0.004 

- 0.021 


Ol = {daSa)v-A{u0Ufj)v-A — {UaSa)v-A{d(3Uf3)v-A 
O 2 = {daSa)v-A{upUi3)v-A + {UaSa)v-A{d/3Ui3)v-A 

+2{daSa)v-A{dpdii)v-A + 2{daSa)v-A{s0Sii)v-A 

O 3 = 2{daSa)v-A{U/3Ui3)v-A + 2{UaSa)v-A{.dl3Up)v-A 

— {daSa)v-A{di3di3)v-A — {daSa)v-A{S/3Si3)v-A 
O 5 = {daSa)v-A{u/3UfS + d/jd/j + Sf^Sp)v+A 
Oa (^daS p')y— j^i^UpUa 'H df^dfy SfjSa^y. 


( 4 ) 


(5) 

( 6 ) 

(7) 

( 8 ) 
(9) 


The values of the coefficients are taken from ref. |^. We choose the version with the flavor 
dependent Aqcd, rrit = 200 GeV/c^ and po = 0.24 GeV. The value of po is chosen so as to 
give as(/ig) = 1 for Aqcd = 0.1 GeV. Among the 4-quark operators above, O 3 contains a part 
which induces the AJ = 3/2 transition, while the others are purely AJ = 1 / 2 . One sees that 
the Ol component is enhanced, which is purely AJ = 1/2. 

It is known that this AJ = 1/2 enhancement alone cannot explain the observed ratio of 
AJ = 1/2 and AJ = 3/2 for the non-leptonic decays of JJ, A and other strange hadrons. 


6 







Several studies have shown that various non-perturbative effects at low energy are crucial in 


understanding the large AJ = 1/2 enhancement |]T^,|T^,|T^ . 


3. Direct quark mechanism 

The weak transition of two baryon systems can be described by a transition potential V. 


We evaluate V in the hrst order perturbation theory in the effective weak Hamiltonian -H/jj ^ 
derived in the previous section, 



( 10 ) 


We employ the quark cluster model for the quark component of the two baryon systems, 


BB'(L, S, J)) = ^®|0(123)0(456)x(L, J)) 


( 11 ) 


where 0 is the internal wave function of the baryon in the non-relativistic quark model and 
x{R) is the wave function for the relative motion with R, the relative coordinate of two 


baryons|p.6||. The operator antisymmetrizes the six quarks. The transition potential eq. (|T^) 


describes the direct quark processes, in which all possible exchanges of quarks between two 
baryons are included (Fig. ^). These direct quark processes are independent from the meson 
exchange diagrams in our formulation because the non-relativistic formalism does not allow a 
pair of quark and antiquark in the intermediate state. Therefore the full transition should be 
given by a superposition of the direct quark and the meson exchange processes. 

In evaluating the transition potential we make the non-relativistic reduction of 
i.e. the Breit-Fermi expansion to hrst order in p/m. The result is given in terms of a set of 
non-relativistic operators listed in Table |l]. The vectors qij and Pi are dehned by 


Qij =^i-Pi= Pj - 



( 12 ) 
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Table 1: The set of transition operators 


Alij — 

[{dh)i{u^u)j + 


1 ® 

1 

A2ij = 

+ 


{di ■ dj) (g) 

1 

A3ij = 

+ 


{di - dj) ® 

{^ij) 

AAij = 

+ 


0 

{P, - P,) 

A5ij = 

— 


0 

{P + P,) 

A6ij = 

+ 

(g)i {di X dj)(g) 

{^ij) 

Alij = 

+ 



{P - P,) 

A3ij = 

— 



(P + P,) 

A9ij = 

— 


{di + dj) (g) 

{d.ij) 

A10,j = 

— 



(P - Ps) 


+ 


0 

{P + g) 

Blij = 

[{d^s)i{d'^d)j + {d^d)i{d'^s)j]® 

1 (g) 

1 

Bllij = 

+ 


{di + dj) (g) 

(P + P,) 


C\ij =[{u'^s)i{Su)j + {d'^u)i{u^s)j]® 1 0 1 


+ 


® {ffi + aj) (g)(Pi + Pj) 



N 


N 


N 


N 


A 

N 


N 

N 








Figure 2: The diagrams for the direct quark mechanism 

where pi denotes the momentum of the i-th quark. In Table the color operator is suppressed, 
that is unity. Among those operators, the operators Al, A2, Bl, B2,C1 and C2 are parity 
conserving, while the others are of hrst order m p/m and parity violating. These operators 
are symmetric in subscripts i and j. Appendix A gives the explicit form of in terms of 

J2i<j ~ Cllij. Because we truncate the expansion at p/m, the change of the orbital 
angular momentum, AL, is restricted to 0 or ±1, namely no tensor transition is allowed. 
In the present study, we restrict our initial state to L = 0 and 1. Table |] shows 24 possible 
combinations of L, S, J, and I for the initial and hnal states. Note that the transition between 
^Pi and ^Pi vanishes because the spin change operator should change the parity as well. 
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Table 2: Possible initial and final quantum numbers for transitions with initial L = 0 or 1 


Ch 

isospin 

spin 

orbital 

Type 

W 

1 

pA pn 

'^0 


1 

1 

2 



^^Po 

2 

-V3 

3 




3 

1 

4 



^ ^Pi 

4 

1 

5 




5 

-V2 

6 


^Pl 


F 

1 

7 



^iPi 

G 

1 

8 



-^3^1 

F 

fb 

V 2 

9 


^Po 


B 

-V3 

10 



^^Po 

A 

1 

11 


3Pi 


C 

-V2 

12 




A 

1 

13 




C 

2 

14 


^P2 

^^P2 

A 

1 

15 




B 

fz 

V 2 

16 



-^3^2 

C 

3 

2 

17 

nA —> nn 



6 

1 

18 



^^Po 

7 

-V3 

19 



^^Pl 

8 

-V2 

20 


^Po 


I 

-V3 

21 



^^Po 

H 

1 

22 


3Pl 

^3p^ 

H 

1 

23 


^P2 

^3p2 

H 

1 

24 




I 

fz 
V 2 
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Matrix elements of ~ are calculated for the quark cluster model 

wave function. Under the condition that |00x) is totally antisymmetric for the exchange 
{l,2,3}<->'{4,5,6}, the matrix element, for instance 




(13) 


is equal to 


^(#xl 411 i2 |#x) + ^(00x1 4 II 36 l#x) 

- ^( 00 x 1 4II12P36 | 00 x) - ^( 00 x 1 4II13P36 | 00 x) 

36 36 

- ^(00X1 41125^36 |00X) - ^(00X1 41135^36 |00X) 

- ^(00X1 41136^36 |00X) 


(14) 


where P 36 represents the permutation operator, 3 <->■ 6 . Fig. ^ shows the diagrams correspond¬ 


ing to each term of eq.(M). Because each baryon wave function is totally antisymmetrized, 
of the initial and hnal state can be replaced by a single P 36 operated to the initial state. N 
is the normalization factor, which depends on the channel but is nearly equal to 1 in general. 
We factorize each matrix element in eq.([T^) into the flavor-spin, orbital and color parts as 


(00x(P, F/, S'/, J)| AU/(P36) |00x(fc, Li, Si, J)) 

{ 00 xflavor-spin| ^^flavor-spin ^^flavor-spin^ j^^^flavor-spin^ 
X (00xorbital| ^iorbital(porbital) |^^^orbital^ 

X (00x^01°’^ I |00x^°^°’^) X W . 


(15) 


W is an algebraic factor required when we factorize the spin and orbital matrix elements. It 
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is defined by 


W = (-)V2A + lV2J + l< 


5,; J 


A° A" 0 


2Lf + l 


^2Sf + 1 


(L„ A«, L„ 0|L^, L,) {Si, A^ 0|5;, 5,) 


(16) 


U/ jj 

where A° and A® are the ranks of the orbital and spin operators respectively, and A = A° = A^. 
In Table we label the combinations of the initial and the final spin-flavor part of |</>0x) 


by “Type”. They are listed in Table |I^ in Appendix B. In evaluating the flavor-spin matrix 
element, we use the SU{6) flavor-spin wave function for the nucleon and A. 

Appendix C is devoted to evaluation of the orbital matrix elements 


where ooi'bital jg following operators. 


(17) 


lij X 5(^. + jfj - Pi- Pj) (18) 

Qij X <^(^i + l^j -Pi- Pj) (19) 

{Pi -Pj) X 6 {]^i + -Pi- Pj) (20) 

{Pi + Pj) X 6 (^i + -Pi-Pj) . (21) 

Color matrix elements are given by 

(color-singlet 111 color-singlet) = 1, (22) 

(color-singlet I |color-singlet) = 1/3. (23) 


4. Induced transition potential 


The obtained AN — 

V{k, k') Li,Si,j 
Lf,Sf,J 


NN transition potential is written in the following form. 

= f(k, kj, + V‘ g{k, kjt + Vj h(k, kj,} X W 

V 2 


(24) 
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The coefficients, and and the functions, /, g and h are given in Tables ^ o 

The numbers are normalized in the unit of 1/660 x v/6/1296 for Type 1, 2, 3, 4, 5, A, B, C 
and F and in the unit of 1/660 x a/3/648 for Type 6, 7, 8, H and I. In Tables ^ ~ Lf), 

exp[ij] and F{L,X) are dehned by 


5{L,,Lj) 
nL^,Lf) 
exp [12] 
exp [13] 
exp [25] 
exp [35] 
exp [36] 
F(0,X) 
F(1,X) 


hiLf 


coseY^.{n)dn 


3^6 

24 
121 
3^3 


exp 

exp 


-b^-{k^ + k'^) 


exp 


—(7fc2 + 13fc'2) 

33^ ^ 


6 


24v^ 


121 


■exp -b^ — {13k^ + 7k'^) 
33 


= exp 


-b^-{k^ + k'^) 

3 


= 47r 


smh{X b'^kk') 


= —Att 


Xb‘^kk' 

f smh{X b'^kk') cosh.{X b‘^kk') 


X%^k‘^k'^ 


Xb^kk' 


(25) 

(26) 

(27) 

(28) 

(29) 

(30) 

(31) 

(32) 

(33) 


The transition potential depends on three quark model parameters that are m, the constituent 
u, d quark mass, m^, the strange quark mass and the Gaussian parameter b. The ratio rris/rn 
is chosen as 5/3 in the calculation of and V^. We use m = 313 MeV and b = 0.5 fm 

in the following section. 

In order to study the contribution of AI = 3/2 component of the we also calculate 

the transition potential without the AI = 3/2 component. Table I Table give the 
coefficients when we omit the AI = 3/2 component. 
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Table 3: The coefficients V-^ for the full Hamiltonian. 


Ty. 

v/ 

h/ 

^3 

v/ 

v/ 

^6 

v/ 

1 

-77588.4 

29011.5 

-787.2 

6117.8 

-10195.3 

39068.5 

-34345.4 

3 

-77588.4 

29011.5 

-787.2 

303.9 

-10195.3 

4185.1 

538.0 

G 

-77588.4 

-79950.0 

-163442.3 

-79950.0 

-50544.8 

-162261.5 

-162261.5 

A 

-77588.4 

-44941.4 

-34614.9 

-13581.0 

-22070.8 

-13756.8 

-12855.1 

6 

-77588.4 

29011.5 

-787.2 

-67.5 

-10195.3 

1956.7 

2766.4 

H 

-77588.4 

-44941.4 

-34614.9 

-13581.0 

-22070.8 

-9633.2 

-16978.7 


Table 4: The coefficients for the full Hamiltonian. 


Ty. 




VI 

VI 

^6^ 

VI 

2 

-54.3 

-297.1 

15321.8 

5760.4 

-43.6 

21392.9 

-2476.8 

4 

-54.3 

-84.6 

69818.0 

31586.4 

1490.8 

85572.9 

84357.2 

5 

54.3 

-122.8 

-14966.7 

-5465.2 

-1767.2 

-758.7 

-2202.1 

F 

-54.3 

184.6 

314.9 

-256.9 

-6683.7 

-57943.0 

-52750.3 

B 

-54.3 

214.6 

166.3 

-2266.7 

-2239.4 

-25126.8 

12162.8 

C 

54.3 

45.2 

-142.0 

53.9 

-182.6 

-1155.6 

-261.3 

7 

-54.3 

115.3 

15321.8 

5760.4 

656.1 

6823.1 

12093.1 

8 

54.3 

152.1 

-14941.7 

-5465.2 

-1467.3 

66.0 

-3026.8 

1 

-54.3 

-60.3 

216.3 

207.4 

-2039.5 

-3959.3 

-9004.7 
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Table 5: The coefficients for the full Hamiltonian. 


Ty 

If 





T.rh 

*^6 


2 

-54.3 

-297.1 

-5116.0 

-4602.0 

-1294.7 

-28090.4 

17845.1 

4 

-54.3 

-84.6 

-23274.8 

-759.8 

-4683.2 

-43698.3 

-41266.8 

5 

54.3 

-122.8 

5007.7 

-586.3 

4817.9 

-10152.2 

-6002.8 

F 

-54.3 

184.6 

0.0 

13521.8 

19700.0 

110914.1 

106694.0 

B 

-54.3 

214.6 

69.6 

1436.8 

6307.2 

21951.0 

1772.9 

C 

54.3 

45.2 

95.6 

295.2 

439.2 

3666.7 

584.1 

7 

-54.3 

115.3 

-5116.0 

552.4 

-1469.7 

-2249.6 

-7995.7 

8 

54.3 

152.1 

4907.7 

444.6 

4742.9 

-8502.8 

-7652.2 

I 

-54.3 

-60.3 

-130.4 

1436.8 

6257.2 

9305.5 

14418.4 


Table 6: The functions /j 

f(k,k'),= i i5(Li,L,)V^^~S(k' -k) 

f(kXh = ~\ 6(Li,Li) exp[12] f 

nk, k% = ~l SiU, Li) exp[13] F (i., 

36 1 

f{k, k'), = -^ 3 S{L,, Lf) exp[25] F (L„ 0) 
f{k, k% = -^l Lf) exp[35] F 
f{k,k\= ^ 1 5{L„Lf) exp[36] F 

f{k, k')r = ^ 5{Li, Lf) exp[36] F (^Li, 0 
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Table 7; The functions gi 


g{k, k')i 
g{k, k ')2 
9{k, k')3 
gik,k')i 
gik, k')^ 

9ik,k')(i 
g{k, k')r 


h{k,k')i 
hik, k')2 
h{k, k% 
h{k,k')4 
h{k,k'), 
h{k,k\ 
h{k,k')j 


^(^4) 

Qfi 1 k 

“iV 3 ~ exp[ 25 ] F(L/, 0 ) 

^ 1/(Li,L/) b expl36| -F 


Table 8: The functions hi 

^(^4) 

3fi 1 k' 

--- I(U,L,) - exp|25] F(L,,0) 
^ 1 I(Li,Ls) h exp[36] -F(ii,|[) 

^ 4 ’ 5 ) 
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Table 9: The coefficients for the AJ = 1/2 Hamiltonian. 


Ty. 

vi 

vi 


v/ 

v/ 

^6 

v/ 

1 

-77588.4 

29011.5 

-787.2 

1994.3 

-10195.3 

14327.3 

-9604.2 

3 

-77588.4 

29011.5 

-787.2 

303.9 

-10195.3 

4185.1 

538.0 

G 

-77588.4 

-79950.0 

-163442.3 

-79950.0 

-50544.8 

-162261.5 

-162261.5 

A 

-77588.4 

-44941.4 

-34614.9 

-13581.0 

-22070.8 

-11007.7 

-15604.2 

6 

-77588.4 

29011.5 

-787.2 

1994.3 

-10195.3 

14327.3 

-9604.2 

H 

-77588.4 

-44941.4 

-34614.9 

-13581.0 

-22070.8 

-11007.7 

-15604.2 


Table 10: The coefficients for the AJ = 1/2 Hamiltonian. 


Ty. 




Vi 

Vi 


Vi 

2 

-54.3 

-22.2 

15321.8 

5760.4 

422.9 

11679.7 

7236.5 

4 

-54.3 

-84.6 

69818.0 

31586.4 

1490.8 

85572.9 

84357.2 

5 

54.3 

60.5 

-14950.0 

-5465.2 

-1567.3 

-208.9 

-2751.9 

F 

-54.3 

184.6 

314.9 

-256.9 

-6683.7 

-57943.0 

-52750.3 

B 

-54.3 

31.4 

199.6 

-617.3 

-2106.2 

-11015.1 

-1948.9 

C 

54.3 

45.2 

-142.0 

53.9 

-182.6 

-1155.6 

-261.3 

7 

-54.3 

-22.2 

15321.8 

5760.4 

422.9 

11679.7 

7236.5 

8 

54.3 

60.5 

-14950.0 

-5465.2 

-1567.3 

-208.9 

-2751.9 

1 

-54.3 

31.4 

199.6 

-617.3 

-2106.2 

-11015.1 

-1948.9 
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Table 11: The coefficients for the AJ = 1/2 Hamiltonian. 


Ty. 




vt 

Vi^ 

T/h 

*^6 


2 

-54.3 

-22.2 

-5116.0 

-1165.8 

-1411.3 

-10863.2 

617.9 

4 

-54.3 

-84.6 

-23274.8 

-759.8 

-4683.2 

-43698.3 

-41266.8 

5 

54.3 

60.5 

4941.0 

101.0 

4767.9 

-9052.6 

-7102.4 

F 

-54.3 

184.6 

0.0 

13521.8 

19700.0 

110914.1 

106694.0 

B 

-54.3 

31.4 

-63.7 

1436.8 

6273.8 

13520.7 

10203.2 

C 

54.3 

45.2 

95.6 

295.2 

439.2 

3666.7 

584.1 

7 

-54.3 

-22.2 

-5116.0 

-1165.8 

-1411.3 

-10863.2 

617.9 

8 

54.3 

60.5 

4941.0 

101.0 

4767.9 

-9052.6 

-7102.4 

I 

-54.3 

31.4 

-63.7 

1436.8 

6273.8 

13520.7 

10203.2 


5. Application to light hypernuclear decays 

We apply the transition potential to non-mesonic weak decays of light hypernuclei. We 
assume that the decay of A in nucleus is incoherent and that one can neglect hnal state 
interactions for two energetic outgoing nucleons and interference effects arising from the an- 
tisymmetrization of the hnal state. Then the decay rate of a light hypernucleus is given by 
a sum of two-body AN NN transition rates. The decay of ^He, for instance, can be 
described in terms of the spin averaged two-body transition rates, TAp^pn and TAn^nn, as 

T ^A^e) = 2 TAp^pn + 2 TAn^nn • (34) 


where 


T = - V 
4 ^ 

Si,Mi 


d^K 


2t^5{E.C.) 


d^p f d?q 


(27r)^ J (27r)' 


i)fin{q]K) 1/Si.Mi (p, ^ 


(35) 
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Here and 'ipfin{q;K) are the initial and final two-body wave functions and V{p,q) is 

the transition potential. In the present study, we also neglect the binding energy of the initial 
state, and use the following energy conservation rule. 


HE.C.) = !^S[K - K- 


(36) 


where K* = 415.9 MeV satisfies 


-\- M]^ — 2 


K* 


2M. 


-\- 2M]sf 


N 


(37) 


Decomposing -0 into partial waves and performing the K integration, one obtains 


r = 


1 MnK* 1 


(27r)2 2 4 


7 E 


L^ ,S,J,m 




(38) 


where 'ipi is the radial part of the wave function and V Li,Si,j (p, q) is the partial wave decom- 

Lf,Sf,J 

posed transition potential. For the decay of S-shell hypernuclei, we neglect L ^ 0 components 
of the ground state wave function and thus only the transitions from the S wave initial states 
to the final plane wave are considered. They correspond to the channels, 1 ~ 5 for Ap —> pn, 
and 17 ~ 19 for An —>■ nn, in Table We label the transition amplitudes by a through / as 


dp d^ 


bp byi 


dn 


fp fn 


1 C _. 1 C 

On —>■ On 




according to the widely used notation |]^. Among them, the amplitudes a, c, and d describe 
the parity conserving transitions, while the others violate parity. By writing the amplitudes 


If = 1 


1 

0 

0 

0 

1 
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simply as 


A'’)Kr'”‘"‘(p. 9)V'r(p) (39) 

for instance, one obtains 

^Ap^pn = 2(27r)^ + l^pl^ + 3|cpp + 3\dp\‘^ + 3|epp + 3|/pp) (40) 

+ 1 ''”''+ 

Note that the 1 = 0 states are allowed only for Ap pn while the J = 1 hnal states are 
allowed both for An —> nn and Ap pn. 

We employ simple wave functions for the initial and hnal states. We use the Gaussian 
with a short-range correlation function for the initial state, 

'ipiniiR) = N^g{R) exp (42) 

where g represents the short range correlation, 

— r R^ 

g{R) = 1 — C exp-^ . (43) 

L ^0. 

For the hnal state we use the plane wave with the same short range correlation function, 

K*) = 9{R) exp {iK* ■ i?} . (44) 

We use the same g{r) for the initial AN and the hnal NN only for simplicity. We choose 

the parameter B as \/2 x 1.3 fm for the S-shell hypernuclei, which corresponds to the shell 

model wave function with the Gaussian parameter 1.3 fm for both the nucleon and A. For 
the short range correlation we choose C = 0.5 and tq = 0.5fm = 6 in the present calculation. 
The strength C is chosen arbitrarily, while we hnd that the results are qualitatively the same 
for other values of G except for the proton asymmetry parameter oi (see section 6). 
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The one-pion exchange (OPE) amplitudes are also computed for the same wave functions. 
We take the form factor p into account, 

M = JJJ d^Rd¥id% - y)^OPe(o - k'2)p(h’2 + , (45) 

where 


p(r) = exp 


3 -Q 

2V 


(46) 


represents the quark density in the baryon. Because AJ = 1/2 is assumed for the weak AAtt 
vertex, OPE amplitudes satisfy an/dp = hn/hp = fn/fp = \/ 2 , while the AJ = 3/2 amplitudes 
give the ratio — 1/a/2. The weak ANti coupling constant is adjusted to the free A decay rate. 


The results for the two-body transition amplitudes are listed in Table |^. The numbers 
given under “AJ = l/2(3/2)” are the results with the pure AJ = l/2(3/2) transition potential. 
We hnd that the direct quark (DQ) amplitudes are in general comparable to the OPE ones, 
especially dp, bp, fp and /„ for DQ are larger than those for OPE. Although dn in the full DQ 
is small due to the cancellation of AJ = 1/2 and AJ = 3/2 amplitudes, it seems accidental 
because by changing the short-range correlation factor, the cancellation may disappear. While 
OPE contains only the AJ = 1/2 component, we hnd large AJ = 3/2 contributions for the 
J = 0 DQ transitions, dp, bp, dn and bn- The AJ = 3/2 contributions for fp and /„ are small. 
The DQ amplitude dp is zero because we neglect the tensor operator by truncating the p/m 
expansion at the order {p/m). On the other hand, OPE has a large dp which comes from the 
tensor part of the one-pion exchange interaction. It is enhanced due to the large momentum 


transfer. 
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Table 12: Calculated transition amplitudes in 10 MeV 



isospin 

spin 

orbital 

full 

Direct Quark 

AI= 1/2 

A/ = 3/2 

OPE 

Ojp 

pA pn 

'Ao 


-78.1 

-23.4 

-54.7 

2.2 

bp 




-53.5 

2.0 

-55.5 

-24.8 

Cp 



-^Ai 

-1.0 

-1.0 

0 

2.2 

dp 




0 

0 

0 

-86.8 

€p 




-23.2 

-23.2 

0 

-43.0 

fp 



-^^Pi 

-55.4 

-53.8 

-1.5 

20.2 

O-n 

nA ^ nn 

'Ao 

^'Ao 

5.5 

-33.1 

38.6 

3.1 

bn 




42.2 

2.9 

39.3 

-35.1 

fn 




-75.1 

-76.2 

1.0 

28.6 


6. Decays of S-shell hypernuclei 


The spin averaged transition rates are decomposed as 


rAp—>pn ^ (^pO T Shpi) , 


rAn—>nn ^ (^^,0 T ST^]^) , 


(47) 

(48) 


where T^rj is the transition rate for the two-body AiV system with angular momentum J, 


■p _ MnK 12 I 1^ |2nj 

(l®pl + l^pl ) 


2(27r)2 

MnK\,_ ,2 


2(27r)2 

■p _ ^nK I2 I |2n^ 

hnO — (I®"-1 + Pn| J 


(|cp| -I- |dp| + \ep\ -I- \fp\ ) 


(49) 

(50) 

(51) 
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Tnl — 


MnK* 2 
I Jn \ 


(52) 


2(27r)2 

Our result for Fj are given in Table [l^. We find that the J = 0 proton-induced transition 
rate, Fpo, is strongly enhanced due to the AI = 3/2 transition. Compared with the OPE 
result, Ppo for DQ is much larger and in fact is dominant while OPE is dominated by the 
tensor transition included in Ppi. This dominance of Ppi in OPE makes the n-p ratio, Rnp, 
small, where 


(53) 


(54) 


■p 

p — ^ neutron induced 

^np — p 

^ proton induced 

For the spin-average hypernuclei, iHe, this ratio is given by 

P _ PnO + 3Pnl 

^np — op • 

i pO -h Oi pi 

In DQ, P„i is also large so that the spin averaged Rnp is as large as 1. Thus we find that DQ 
and OPE predict qualitatively different values for i?„p(^He), while the decay rates, P(^He), are 
roughly equal. The experimental data prefers DQ, which indicates a signihcant contribution 
of Pnl. 

Recently, Schumacher proposed to check the AI = 1/2 rule in the non-mesonic decays of 
the S-shell hypernuclei |T^ , |2D| . He calculated the ratios of P nj by using the following relations 
and the corresponding experimental data, 

PnO + 3P„i 


RnpilHe) = 
RnpiiHe) = 


Ppo -|- 3Ppi 

2rn0 


Ppo + 3Ppi 
P(^He) _ Ppo -|- 3Ppi -|- 2P„o 


(55) 

(56) 

(57) 


^(aH) 2PpO -|- PnO + 3Pnl 
where it is assumed that the Patj’s are common for ^He, ^He and ^H. These equations 
determine three ratios of Patj’s. Using the experimental values of PAp^pn and PAn^nn, we 
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Table 13: Decay rates of light hypernuclei. All the decay rates are in the unit of T free- The 
experimental data for TAp^pn, ^An—nn.) r^He) and i?„p(AHe) are taken from ref. [18]. Those 
for RnpiXRe), T^ ,„.(AHe)/r„ ,„.(aH) and T^o/Tpo are taken from ref. [19][20]. See the main 
text text for the “ experimental” T atj’s. 



Direct Quark 

OPE 

DQ ± OPE 

Exp 


Full 

A/=i 


+ 

- 


Tpo 

0.177 

0.010 

0.012 

0.235 

0.143 

0 ~ 0.116 

Tpi 

0.071 

0.067 

0.193 

0.260 

0.269 

0.074 ~ 0.187 

TnO 

0.035 

0.021 

0.024 

0.002 

0.118 

0.063 ~ 0.553 

Tnl 

0.111 

0.114 

0.016 

0.042 

0.212 

0.049 ~ 0.196 

r Ap—^pn 

0.097 

0.053 

0.143 

0.253 

0.238 

0.105±0.035 

r An—>nn 

0.092 

0.091 

0.018 

0.032 

0.189 

0.100±0.055 

raHe) 

0.378 

0.295 

0.333 

0.573 

0.854 

0.41 ± 0.14 

ii’„p(^He) 

0.94 

1.70 

0.12 

0.12 

0.79 

0.93 ± 0.55 

Rnpiiile) 

0.18 

0.20 

0.08 

0.004 

0.24 

0.18 ± 0.12 

r„.^.(4He) 

0.63 

0.66 

6.58 

1.69 

1.14 

1.65 ± 0.77 

Tno/rpO 

0.20 

2.00 

2.00 

0.01 

0.854 

-0.8 ±2.7 

ai(iHe) 

0.01 

0.02 

-0.19 

0.20 

-0.44 
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evaluate the “experimental” values of T^j given in Table O The ratio T^o/Tpo is especially 
sensitive to the AJ = 3/2 mixing, he., it is 2 for the pure AJ = 1/2 transition, while it 
becomes 1/2 for the pure AJ = 3/2 transition. The second row from the bottom in Table ^ 
gives the ratio T^o/Tpo- DQ gives a much smaller value than 2, which clearly demonstrates 
the contribution of AJ = 3/2. The present data, —0.8 ± 2.7 for the S-shell hypernuclei are 
not conclusive. One also sees that DQ mechanism can reproduce the n-p ratio for ^He, and 
the ratio r(^He)/r(^H) fairly well. 

When the hypernucleus is polarized, the angular distribution of the outgoing proton has 


an asymmetry. It is parameterized in terms of the asymmetry parameter oi as 


W{0) = l + ai(p)PAA(cos0) . 


(58) 


where Pa is the polarization of A and 6 is the angle of the outgoing proton to the A polarization. 
This parameter for |He is given by||TT 

2 ^ + dp)fp 


ai(iHe) = 


(59) 


al + bl + 3(c2 + dl + el + f^) ' 

Recent experimental data indicate a large negative ai(p) for p-shell hypernuclei, ai(p) < 
— 0.6|^. Our result for ai(|He) is very small because our dp is zero and Cp is also small. But 
the result is rather sensitive to the choice of the short-range correlation, and therefore is not 
conclusive. 

So far we have not considered the interference of the DQ and OPE amplitudes. As we 
have argued in section 3, the present formalism allows us to regard OPE independent from 
DQ and therefore to superpose these two amplitudes. Because the relation between the 
phenomenological ANtt vertex in OPE and the effective weak Hamiltonian JJ/jj=^ in DQ is 
not known, the relative phase of the two amplitudes cannot be determined. Thus we evaluate 
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DQ ± OPE and the results are listed in Table O. One finds that the difference between the 
two choices of the relative phase mostly appear in the neutron-induced decay rates. Pnj’s are 
suppressed in (DQ -|- OPE) and thus the n-p ratio Rnp becomes very small. In this sense, the 
experimental data prefer the (DQ — OPE) combination. (DQ — OPE) also predicts a large 
negative ai(^He), which seems to agree with the experimental value for the p-shell hypernuclei. 
The ratio Pno/Ppo tends to be small (-C 2) for both ( DQ ± OPE ) and again indicates a 
large AJ = 3/2 contribution. In both (DQ ± OPE), we find that PAp^pn. is overestimated 
and therefore the total decay rate P(^He) is too large. This is again due to the large tensor 
component in Ppi(OPE). It seems important that OPE is calculated in the quark interaction 
point of view in order to make a reliable prediction for the DQ - OPE interference. 

Recently Ramos and Bennhold studied contribution of the heavy meson exchanges such 
as K, p and rj . They indicate that such contributions are suppressed by the short range 
correlation and furthermore, they tend to cancel with each other. Recent studies of the 27r 
exchange mechanism indicate that the diagrams with SA and NN intermediate states cancel 
with each other and the net effect contributes only to the J = 0 amplitudes. In all, the 
meson exchange contributions other than OPE seem to be small. Therefore one may describe 
the AN —NN transition well only by the DQ and OPE. 


7. Discussions and Conclnsion 

Nonmesonic decays of hypernuclei provide us with a new type of the hadronic weak in¬ 
teraction. The large momentum transfer (due to the mass difference of A and N) makes the 
transition sensitive to the short distance quark structure of the two baryon system. Indeed, 
it is found that the contribution of the direct quark processes is as large as that of the con- 
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ventional one pion exchange weak interaction. Furthermore, we have found that the J = 0 
transition amplitudes show a large AJ = 3/2 contribution and therefore that the AJ = 1/2 
rule is signihcantly broken. This may be the hrst clear evidence for the AJ = 3/2 weak 
transition, that is expected in the standard theory of the weak interaction. 

We have employed, in the present analysis, an effective weak Hamiltonian for quarks, 
which takes account of the one-loop perturbative QCD corrections. Then we have evaluated 
the transition amplitudes using the quark model wave functions of baryons to the hrst order in 
the weak interaction. The havor/spin structure of the amplitudes rehects the SU(6) symmetry 
of the baryon wave functions, which have been verihed in the low energy baryon spectrum 
and properties of the baryons. 

We have derived an effective AN NN transition potential, and applied it to the s-shell 
hypernuclear decays. It is found that the decay amplitudes show distinctive features when 
they are compared to the one-pion exchange. Especially, the ratio of the neutron-induced 
and the proton-induced decay rates is discriminative of these mechanisms. It is suggested 
that the ratios of the transition rates with various spin-isospin specihcation can be obtained 
from the experimental data for the s-shell hypernuclei and they are useful in testing different 
mechanisms of the transition. Further experimental studies are most desirable. 

There are a number of remaining problems. The relation between the phenomenological 
A —>■ Ntt transition Hamiltonian and the effective quark Hamiltonian is to be studied. It is 
favorable to apply the same quark Hamiltonian to the mesonic decay as well so that a unihed 
view of the hypernuclear decay is obtained. The AJ = 1/2 enhancement mechanism for the 
mesonic decay is especially important in this regard. This line of study is underway and will 


be reported elsewhere[25 . 
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For hypernuclei other than the s-shell systems, we need a realistic calculation combined 


with the nuclear structure analysis. We have provided the baryonic two-body transition 
potential that can be used in any hypernuclear structure calculations. Because of the nonlocal 
structure due to the quark exchange effects the transition potential is given in the momentum 
space, but the transformation into the coordinate space is straightforward. 

We have not considered so far the second order process with a S — intermediate state 
induced by a strong pion (meson) and/or quark exchanges. The weak EiV —>■ NN decay can 
be also computed in the same direct quark mechanism. It is found that the mixing of T^N does 
not change our main conclusions mentioned above, though its contribution is not negligible 


quantitatively. The results of this calculation will be published in a separate article pq . 





Appendix A. Non-relativistic forms of 0i ~ O, 


The Hamiltonian ^ is given by eq.(|D: 

(a ~ fio) =-% E km (60) 

The operators, O 2 ~ Oq, contain the terms {ds){ss), which we omit because they do not 
contribute in the valence quark model. Then the Breit-Fermi expansion of Oi ~ Oq is given 
in terms of the operators Al ~ Cll by 

Oi = (A1 - Cl) - (A2 - C2) - 2(^3 - C3) 

- C4) - 1{A5 - C5) 

-‘^^(A6 - C6) - 1(A7 - Cl) ~ 1(A8 - C8) (61) 

O 2 = {A1 + 2B1 + Cl) - (A2 + 2B2 + C2) - -(A3 + 2B3 + C3) 

+ 254 + CA) - ^(A5 + 255 + Cb) 

- 2/i + ^ (A 6 + 256 + Cb) -^iA7 + 257 + C7) - ^( 2 l 8 + 258 + C8) (62) 

O 3 = (2A1 - 51 + 2C'l) - (2A2 - 52 + 2C2) - ^(2A3 - 53 + 2C3) 

- 2/i + <^ ( 2 A 4 _ 54 + 2(74) - ^(2A5 - 55 + 2C5) 

- ^^^^ (2A6 - 56 + 2C6) - ^(2A7 - 57 + 2C7) - ^(2A8 - 58 + 2C8) (63) 

err 

O 5 = (A 1 + 51) + (A 2 + 52) " 2(^6 + 56) + -(A7 + 57) + -(A 8 + 58) 

~{A9 + 59) - ^^^^(AlO + 510) + ^(All + 511) (64) 

{ r r r 

(51 + Cl) - -(53 + C3) + -(54 + C4) + -(55 + C5) 

+ ^^^(59 + C9) + ^(510 + CIO) + ^(511 + Cll)| (65) 
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where 


rUs + m , ^ nis — m 

/i = —- and h = —- , 

Zrusm ImsTn 

are given in terms of m, the mass of constitnent u, d qnarks, and rUg, the mass of constitnent 
s qnark. In the present calcnlation, we use m = 313 MeV and m/irtg = 3/5. Note that we 
have made the following Fierz transformation on Oq so that the color part becomes unity. 



Oq {,dQi^Sp)y—j\{UijUQ^ df^dfgi "h 


(67) 


= —‘^{daUa)s+p{upSfs)s-P — 2{dada) S+P^dfjS ff) S-P 


( 68 ) 


Appendix B. Flavor-spin part of |</>0x) 


Table ^ gives the flavor-spin part of |</>0x) for each Type. One sees that |00x) is totally 
antisymmetric under the exchange {l,2,3}'«-*-{4,5,6}. We choose either = 0 or 1 so that 
(sj, A^, S 2 , 0|s/, is not zero. These “Type”s are refered to by the channels in Table ^j. 


Appendix C. Orbital matrix elements 


The orbital part of the internal wave function of A and N is taken as 


^(l_2,3)Orbital _ 


(69) 


where b denotes the size of the baryon. We choose b = 0.5 fm. It is convenient to use the 


Jacobi coordinates, dehned in Table in writing the orbital wave function of the two baryon 


system, 


0(1, 2, 5, ('n+|+fa _ r4+|±^^ 


(70) 
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Table 14: Flavor-Spin part of |00x)- 



Initial 

Final 

Ty. 

Orb. 

Flavor 

Spin 

Orb. 

Flavor 

Spin 

1 

S 

■^{pK + Kp) 

7^(n-iT) 

S 

-^{pn + np) 

^(TT - TT) 

2 


+ 

- 

P 

+ 

+ 

3 


- 

+ 

s 

- 

+ 

4 


- 

+ 

p 

- 

- 

5 


- 

TT 

p 

+ 

TT 

A 

P 

+ 

TT 

p 

+ 

TT 

B 


+ 

^(Tl + IT) 

s 

+ 

Tfd'T “ ■I'!') 

C 


+ 

TT 

s 

- 

TT 

F 


- 

TjdT “ TT) 

s 

- 

■^(TT + TT) 

G 


- 

- 

p 

- 

- 

6 

s 

^(nA-h An) 

^(TT - TT) 

s 

nn 

^(TT - TT) 

7 


+ 

- 

p 

nn 

+ 

8 


- 

TT 

p 

nn 

TT 

H 

p 

+ 

TT 

p 

nn 

TT 

I 


+ 

^(TT + TT) 

s 

nn 

^(TT - TT) 
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Table 15: Definition of the Jacobi coordinates and their conjugate momenta 


62 

= n - 7*2 


P 12 

— j 0 l-jD 2 

2 


■^12-3 

= r3 


P 12 - 

-3 = Ip3 - 

P1+P2 

3 

■^45 

= - r 5 


P45 

_ P4-P5 

2 


■^45-6 

= rg - ^ 


P45- 

-6 = |P6 - 

P4+P5 

3 

R 

_ n+ri+rs rt+Ti+rfi 


P 

_ 'Pl+P2+P3-P4-'P5-m 

3 3 



2 

Rg 

_ ri+r2+D+75+ni+7% 

6 


Pg 

= Pi +P2+P3+P4 + P5+Pe 


The function 0(1, 2, 3 jorbital jg terms of the Jacobi coordinates ^12 and ^ 12-3 as 

(71) 






or in the momentum space as 


0 ( 1 , 2 , 3 )orbital ^ ('g^^2^4/3('g^^2^3/4 {- 52 ^ 2 } 






(72) 


Explicit forms of the orbital matrix elements 


^^(2Tr)H(P' - g')| Oy(Ps6) |#(2ir)“i(P - k) ) 


(73) 


are listed in Table 16 ~ Table 19. The matrix elements 


i/)l Oij ^(Px) \4><j>x(k, L.) ) 


(74) 


are given by the integration 


Id^ dk y^°;(fc') ( 00(27r)3j(P' - p)| |00(27r)3j(P - k) ) Y^k) . (75) 
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Table 16: Orbital matrix elements. 


O^J{P^6) 


{ (p(j){2TTf5{P' 


fc0|O*j(^36)IM2vr)35(P-fc) ) 


ll2 

112 P 36 

113 P 36 

I 25 P 36 

1 35 P 36 

1 36 

I 36 P 36 


sVg 

^exp 

24V^ 


- k) 

+ 12k ■ k' + ISA:'^) 


3^/3 


exp 


-b^^{k^ + k'^) 


24V^ 

121 


■ exp 


exp 

exp 


-b^^{13k^ + 12k-P + 7k'‘^) 

00 

b^hk‘^ + 2k-P + k'^) 

O 

b'^Uk'^ + 2k ■ P + k'^) 

O 


Table 17: Orbital matrix elements. 


O^J{P36) 


{ (^(k{2TTf5{P' 


P)|0,,(P36)IM2vr)35(P-fc) ) 


912 

912 T36 

913 T36 
925 T36 
q3b P 36 
936 

936 T36 


0 

0 

^{-36fc + 12p}(|li3P36|) 
— 33fc + 33fc'| (I I 25 -P 36 I) 
^{-I2fc + 36P}(|135P36|) 
^ |-66fc + 66fe'| (I I 36 I) 
^{ + 22fc + 22p} (I 136^36 1) 
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Table 18: Orbital matrix elements. 


OijiPse) 

{ <l><l>{2nf5{P' - k')\O^.iPse) \cj)(P{2n)H{P - k) ) 

(A - A) 

0 

(A - A) t-36 

0 

(A - A) PsG 

^{+36fe-12fe'}(|li3P36|) 

(A - A) ^36 

— |+33fe + 33fe^| ( I 25 T 36 ) 

(A - A) ^36 

l{_12fe + 36fe'}(|l35P36|) 

(A - A) 

^{+22fe + 22fe'}(|l36|) 

(A - A) P 36 

^ |- 66 fe + 66 fe'| (1 136^36 1 ) 


Table 19: Orbital matrix elements. 


o^,{Pm) 

( <l><l>{2n)H{P' - k')\ 0*,(P36) |</></>(27r)3,5(P - k) ) 

(A + A) 

^{+22k + 22k']{\li2\) 

(A + A) T 36 

(A + A) ^*36 

^{+33fe-33fe'}(|li2P36|) 

^{ + m-48fe'}(|li3P36|) 

(A + A) -P 36 

(A + A) T 36 

0 

^{-48fe + m'}(|l35P36|) 

(A + A) 

0 

(A + A) T 36 

0 
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